QUANTUM CONVOLUTION AND QUANTUM CORRELATION 
ALGORITHMS ARE PHYSICALLY IMPOSSIBLE 



CHRIS LOMONT 

Abstract. The key step in classical convolution and correlation algorithms, 
the componentwise multiplication of vectors after initial Fourier Transforms, 
is shown to be physically impossible to do on quantum states. Then this is 
used to show that computing the convolution or correlation of quantum state 
coefficients violates quantum mechanics, making convolution and correlation 
of quantum coefficients physically impossible. 



1. Introduction 

The Fast Fourier Transform (FFT) is arguably the most important algorithm 
in computer scienee. Many applications, from image processing, signal processing, 
pattern matching, polynomial multiplication, number multiplication, and many 
others are accomplished efficiently by utilizing the FFT to compute a Discrete 
Fourier Transform (DFT) of some set of data (see Chapter 32] and Chap- 
ter 3]). Naively, the DFT of iV = 2" complex values has complexity 0(7V^), but 
the famous paper by Cooley and Tukey |^ introduced the Fast Fourier Transform, 
reducing the complexity to (^{NlogN) = 0(n2"), making Fourier transforms ex- 
tremely useful in computer algorithms. This efficiency is the basis of fast convo- 
lution and correlation algorithms, which do a DFT on each input sequence, then 
a componentwise multiplication, then an inverse DFT, computing convolution or 
correlation with complexity O(iVlogiV). 

In the last decade, quantum computing has become well known due to the integer 
factoring algorithm of Shor ||] and the database search of Grover , both of which 
have complexities much better than their classical counterparts. There are other 
problems [p^ for which quantum computers perform exponentially better than any 
classical (Turing) computer. A quantum Fourier Transform (QFT) can be done on a 
quantum state consisting of = 2" complex values with complexity 0(log^ N) = 
O(n^) (see for example, H, Chapter 5]), making it exponentially faster than the 
classical counterpart, but unfortunately the peculiarities of quantum mechanics 
disallows using this algorithm as a direct replacement for all FFT algorithms. In 
fact; after the breakthrough algorithms of Shor and Grover, no algorithms of similar 
importance have been found, although there has been intense work in this direction. 

Since the QFT and inverse QFT are more efficient than their classical coun- 
terparts, and the FFT and inverse FFT are the cornerstones for convolution and 
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correlation algorithms, it is reasonable to attempt to construct quantum analogues 
of convolution and correlation algorithms that outperform their classical counter- 
parts. The key point of this paper is that there is no physically realizable method to 
compute the normalized convolution or correlation of the coefficients of two quantum 
states. Thus replacing classical algorithms that rely on convolution or correlation 
cannot be done in a simple, direct algorithm replacement manner, but must be ap- 
proached by more sophisticated techniques, avoiding the method examined in this 
paper. The precise definitions of quantum convolution and quantum correlation 
are in section ^, and it is shown in section ^ that they cannot be computed without 
violating quantum mechanics. 

2. Background 

Since this paper touches upon several areas, we review some key ideas to make it 
clearer for readers coming from different backgrounds. In particular, we will recall 
how convolution and correlation are done classically, followed by a review of rules 
for quantum computation, sufficiently covered for our needs. Then we will state 
the convolution and correlation problem for quantum states. Finally, we will prove 
such an approach is doomed to fail, since it would allow a violation of quantum 
mechanics. 

2.1. Classical convolution and correlation. Convolution and correlation have 
many uses in algorithms, which we will not go into here. Since different authors and 
fields of study use slightly different definitions for the following terms, we define 
them here, and choose a form amenable to making quantum versions. In particular, 
different authors use different factors where the V^V appears in our definitions, 
depending on how they will use them. We choose the conve ntions below to make 



our definitions mimic the quantum versions in section 2.2.1 , where the factors are 
dictated by unitarity. 

Definition 1 (Discrete Fourier Transform). Given a sequence ao, ai, . . . , aN-i of 
N complex numbers, the Discrete Fourier Transform (DFT) of the sequence is 
defined to be the sequence ag, ai, . . . , aj^^i given by 

= ^ E "'^-'"^'^'^ ,- = o,i,...,iv-i 

Definition 2 (Inverse DFT). Given a sequence ao, ai, . . . , a^-i of N complex 
numbers, the Inverse Discrete Fourier Transform (IDFT) of the sequence is 
defined to be the sequence ao, ai, . . . , a^v-i given by 

N-l 

i 

3=0 

The IDFT and DFT are inverses, so are 1-1 on sequences. 

Definition 3 (Convolution). The convolution of two sequences of N complex 
numbers Si ~ {aQ,ai, . . .aN-i) and S2 ~ {Po, Pi, . . . is defined to be the 

sequence S3 = (70, Ji, ■ ■ ■ , Jn-i) given by 

N~l 

7fc=^aj/3fc_, for fc = 0, 1, . . . iV - 1 

3=0 



N-l 

^ y a.e-^^*-'"'^/^ fc = 0,l,...,7V-l 
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where subscripts are taken (mod N). 

Definition 4 (Correlation). The correlation of two sequences of N complex num- 
bers Si ~ (ao, ai, . . . aN-i) and S2 = (Po, Pi, ■ ■ ■ Pn-i) is defined to be the sequence 
S3 = (70,71, ■ • ■ ,lN~i) given by 

7fc=^a*/3fe+j for fc = 0, 1, . . . TV - 1 

j=o 

where a* is the complex conjugate of aj and subscripts are taken (mod N) . 

2.1.1. Convolution/correlation algorithm. It is a simple exercise to check that an 
algorithm to compute the convolution (or correlation) of two length N sequences 
in time O(iVlogA^) is: 

(1) FFT - apply the FFT to each sequence in time 0{N log N), giving two 
sequences of length N. 

(2) Multiply - multiply the resulting sequences componentwise in time 0{N), 
assuming each entry can be multiplied in constant time. If doing correla- 
tion, conjugate the first FFT sequence in time 0{N). 

(3) Invert - Apply the inverse FFT to the resulting sequence with time complex- 
ity 0(iVlog A^), and multiply each number by giving the convolution 
sequence. 

2.2. The rules of quantum computing. A quantum system is determined by a 
quantum state vector, often denoted {ip). Mathematically, a quantum state on n 
two-state particles (each called a qubit) is a unit vector in a A^ = 2" dimensional 
complex Hilbert space. For the rest of this paper A^ will be 2" for some positive 
integer n. The Hilbert space is equipped with an orthonormal basis labelled \i), for 
i = 0,1, . . . , N — I. This state is written as 

1=0 

where the Oi are complex numbers, \i) is the i*'* orthonormal basis element, and 
the Oi satisfy the normalization condition joip = 1. Quantum mechanics allows 
two methods to change the state of a system: 

(1) Unitary transformations. Any state change of an isolated system must 
be reversible, and must satisfy the normalization condition above, which 
leads to unitary operations on the state. That is, a state can be 
transformed to the state U\tjj) via the unitary matrix U. Unitary means 
WU = I = UU^ where / is the identity matrix, and f is the conjugate 
transpose. 

(2) Measurement operators. Applying a measurement to a quantum state 
returns an "answer" with a probability related to the magnitude of the coef- 
ficients, and places the system into the state whose value was returned from 
the measurement. Precisely, a quantum measurement consists of a collec- 
tion {Mm} of measurement operators, which must satisfy M^^Mj^ = I 
to preserve probability. When applied to a system with probability 

p{m) = {i;\Ml,AUi,) 



4 



CHRIS LOMONT 



the state becomes the state 




y/pirn) 

and outcome m is observed by the measurement. See [|[ Chapter 2] for 
details and examples. 
(3) Combining states The concatenation of two quantum systems on n and 
m qubits with states {ip) and \4>) respectively is the tensor product state 
(8) |(/)) in complex 2"+™ dimensional space. Shorthand is 



For the purpose of this paper it is necessary to define what we mean by "physi- 
cally realizable computation" . There are many approaches to this question, but for 
our purposes it is enough to take a very general definition, under which we show 
convolution and correlation to be impossible, since they remain impossible under 
more confining and precise definitions. A quantum computation can utilize extra 
states called "ancillary qubits" as working space, so we allow a finite number of 
extra qubits to be used internally to the system. 

Definition 5 (Physically realizable). A mapping P of quantum states 

A 1^) 

is called physically realizable if there exists a finite fixed sequence of unitary 
transformations and measurement operators and a quantum state \p) (called the 
ancillary qubits) performing the mapping. That is, for all |(/)) , the fixed sequence S 
performs the mapping 

\<P)®\p) 

Note in particular, the mapping is a composition of linear operators, so is a lin- 
ear operator. There are some researchers working on nonlinear quantum mechanics, 
but there is no experimental evidence to date that supports nonlinearity One 
interesting result is that if quantum mechanics is nonlinear, then quantum comput- 
ers can solve NP-complete and problems in polynomial time which some 
believe is very good evidence that quantum mechanics is strictly linear. For our 
purposes we use the standard view that quantum mechanics has linear evolution. 

2.2.1. The Quantum Fourier Transform. Similar to the classical Fourier Transform 



and inverse defined in section 2.1, we can define quantum analogues, operating on 



sequences of complex numbers stored as quantum state coefficients. 

Definition 6 (Quantum Fourier Transform). The Quantum Fourier Transform 

( QFT) is the unitary map defined on basis states \j) as 

QFT J_^\2.^3k/N, 



^ fro 



and extended by linearity. 
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Definition 7 (Inverse QFT). The Inverse Quantum Fourier Transform (IQFT) 
is the unitary map defined on basis states \k) as 

\k) ^ _l=^e---^-/^b-) 

and extended by linearity. 



Note 8. It is easy to check these are inverses on quantum states, and that each 
transform is unitary. 

The quantum state \^p) ~ X^'^dO then transforms as 

N-l N-1 
j=0 J 




(2) = >:i^>>,»""""ii*^> 



One reason some quantum algorithms are superior to classical ones is that the 
QFT has complexity O(log'^iV) = 0{n^) j^. Chapter 5], which is exponentially 
faster than the classical 0{n2") version! It is precisely this exponential speedup 
that allows Shor's integer factoring algorithm to out-perform the best known clas- 
sical one. 

3. Problem Statement 
In this section we describe an attempt to construct quantum versions of the 



convolution and correlation algorithms, based on the classical algorithms 2.1.1 
Since the QFT operates on quantum coefficients, we will encode the sequences in a 
quantum state, with the sequence entries stored in the coefficients. Since a quantum 
state ai\i) requires J2i ~ 1j require each sequence to have at least one 
nonzero entry, and then we normalize the sequence to have norm 1. We ignore 
the physical method the quantum state is constructed from the sequence, since the 
many ways to construct it would lead us. So we have a (hypothetical) method of 
converting a nonzero sequence of complex numbers into a representative quantum 
state. 

The problems become: 

Definition 9 (The Quantum Convolution Problem). Given quantum states repre- 
senting the two initial sequences, compute a quantum state representing the convo- 
lution sequence. That is, given the two states (for N = 2" j 

(3) |a) = E"'l*) 

i=0 

(4) I/?) = 

3=0 

compute a state 

N-l 

(5) w^E^'^i^) 

fc=0 
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where \j) represents the normalization of the sequence given by 

(6) Ck=Y.a,(3k-j fork = Q,l,...N-l 

j=o 

and subscripts are taken (mod A^). 

Definition 10 (The Quantum Correlation Problem). Given quantum states rep- 
resenting the two initial sequences, compute a quantum state representing the cor- 
relation sequence. That is, given the two states (for N ~ 2" j 

N-l 

(7) |a) = 5]a,|z) 

(8) m = Y^m 

compute a state 

(9) w^E^'^i^) 

fe=0 

where I7) represents the normalization of the sequence given by 

N-l 

(10) cfe = ^ a*l3k+j fork = 0,l,...N-l 

3=0 

and subscripts are taken (mod A^). 

Note 11. For each definition to make sense, the resulting sequence Cfc cannot be all 
zeros. This can be guaranteed by requiring that the initial sequences are each not all 
zero, and then padding each initial sequence by appending TV zeros. Then, choose io 
and jo each minimal so that 7^ and f3jg 7^ 0. Then Cip+jp — Oii^Pj^ 7^ in the 
convolution sequence. Pick ii maximal so that ai-^ ^ 0, then Cj^^i-^ = a*_^f3jg 7^ 
in the correlation sequence, with subscripts taken (mod N). Padding by N zeros 
requires adding only one more qubit, so we assume for the rest of this paper that 
this is the case. Then each problem is well-defined. This changes each problem 
slightly, but proving impossibility in these slightly amended case suffices to prove 
impossibility in general. 

4. Impossibility proof 

This section gives a proof that the problems given in definitions || and |l^ cannot 
be computed by any device obeying quantum mechanics. Precisely, we prove in 



section 4.3 



Theorem 12 (Impossibility of quantum convolution). There is no physically real- 
izable process P to compute the (normalized) convolution of the coefficients of two 
quantum states. That is, for arbitrary quantum states X^i'^il*) '^'^'^ X]j ^ib); there 
is no physically realizable process P to compute the state 

N~l JV-lJV-l 

(11) ^ a,bj\ii) A ^ ^ ajbk-]\k) 



QUANTUM CONVOLUTION IMPOSSIBILITY 



7 



where A = 1 / v^CI'^^^'TF normalization factor, subscripts are taken (mod TV), 

and N ~ 2" for some integer n > 0. 

and 

Theorem 13 (Impossibility of quantum correlation). There is no physically real- 
izable process P to compute the (normalized) correlation of the coefficients of two 
quantum states. That is, for arbitrary quantum states '^^ai\i) and there 
is no physically realizable process P to compute the state 

N-l N-lN-l 

(12) ^ ^ E E <'*Mi\k) 

where X ~ 1/ y^Yl Wi^jl^ "is the normalization factor, subscripts are taken (mod N), 
and N = 2^ for some integer n > 0. 

This is done in several parts. First, we show that step ^ in algorithm p.l.l| has 
no quantum analogue by studying requirements on a linear transformation that 
attempts step ^. Then we show no physical process consisting of arbitrary sequences 
of unitary operations and measurements can perform step This proves that step 
^ is impossible to compute on quantum states. Finally, we show that any physical 
process able to compute quantum convolution or correlation as in definitions |^ and 
p!o| would be able to compute the impossible step |[ a contradiction. Thus there 
can be no quantum convolution or correlation done on quantum states. 



4.1. Linearity considerations. Suppose we try to obtain the componentwise 
product of two states, using some linear operator on the initial state. Heuristi- 
cally, this should fail, since all components in the outcome may be zero. Even if 
we remove this case (which is reasonable, see note |l^), we show such an operation 
is still not possible. Since linear operators encompass both unitary operations and 
measurement systems, we analyze the linear operator case first. As in definition ^ 
we allow the algorithm to use arbitrary finite "workspace" in the form of a third 
quantum state Ck\k) on m qubits in order to show the problem cannot even be 
computed with extra working space. 

Lemma 14. There is no linear operator L and a quantum state ^ Cfej/c) such that 
for arbitrary quantum states ^ai|i) and X^^jli) with ^ \aihi\ ^ 0, the following 
operation is performed on quantum states: 



(13) ^ ai6jCfc|ijfc) i \Yaihi\i) J ® ^djkljk) 

i,j,k \ i / \ j.k 



where the djk are functions of the ai,bj, and Ck, with X^Mj/cP ~ '^^'^ ^ ~ 
1 / VX); I'^^i^iP ^■s normalization factor. 

Proof. Assume there is a linear operator L and initial state '^k\k) satisfying ^ 
Write L as 

L = ^ erstabc\rst){abc\ 

a,b,c 
r,s,t 
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and apply it to the initial state to get 

( \( 

^^erstabc\rst){abc\ \^^aihjCk\ijk) \ = aihjCkersUjk\rst) 



(14) 



a.h.c 



We want this state to equal the desired outcome, giving 
(15) '^aibjCkerstijk\rst) = A arbrdst\rst) 



r.s,t 



Looking at the \Qst) component, we require for a fixed s, t that 
(16) Xaobodst = ^ aibjCkCostijk 

Equation |l^ must hold for any and bj coefficients giving a quantum state. Fix 
e e (0, 1), and pick sequences of coefficients as 



(17) {aj = {e,yT^,0,0,0...,0} 

(18) {&,} = {1,0,0,...,0} 



Then A = 1/e, ao^o = oi^o = Vl ^ e^, and aifej = for all other i,i combinations. 
Equation becomes 



(19) dst = e^CfeCostoofc + VT~~e2^Cfceostoifc 



(20) = eCi + - £2^2 

where Ci and C2 are constants depending on the initial choices for L and Cfc. 
Normalization requires 

(21) 1 = ^14, 



2 



(22) = ^leCi + yr^CsP 

(23) = MN [e^Cl + 26^1 - e^CiCa + (1 - e2)C|) 

(24) = MNCl + e^MN{Cl-Cl) + 2ey/]~~^MNCiC2 

Equation |2^ must hold true for all e G [0, 1]. When e = 0, the last two terms are 
zero, forcing C2 = l/VMN ^ 0. When e = 1, the last term is zero, and the middle 
term must be zero, forcing C| = C? 7^ 0. Now the middle term is zero for all e. 
However, when e = 1/2, equation |2j is a contradiction, so there can be no linear 
operator L and initial state ^ Cfc|fc) satisfying |l3| for arbitrary input states. 

□ 

Note lemma ^ is a purely mathematical statement, driven by physical require- 
ments, and cast in quantum mechanics language. The following variation is used 
later to prove the impossibility of quantum correlation. 
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Corollary 15. There is no linear operator L and a quantum state ^ Ck\k) such that 
for arbitrary quantum states ^Oiji) and '^bj\j) with ^ \aibi\ ^ 0, the following 
operation is performed on quantum states: 

(25) ^ aibjCk\ijk) — 




with ^ |rfjfeP ~ 1; o,nd X = 1/-^^^ Ifli^ip is a normalization factor. 

Proof. The proof is the same as lemma |l^, since the sequences needed in the proof 
were real numbers. □ 

4.2. Unitary transforms and measurements fail. 

Theorem 16. There is no physical process P capable of either of the transforma- 
tions 

(26) ^aj6j|ij) A^a,;6i|i) 

(27) E^^^^l*^') ^ A^a*6,K) 

on arbitrary quantum states with ^ \aibi\ ^ 0, where X is a normalization factor. 

Proof. Since unitary transforms are linear operators, both lemma |l^ and corol- 
lary |l^ remain true when the phrase "linear operator L" is replaced with "unitary 
transform L" . So no single unitary operation L can perform either transformation. 
Similarly, a measurement operator {Mm} is a collection of linear transforms on a 
quantum state, so no single measurement operator can perform either transforma- 
tion. Finally, given any (finite) sequence of unitary operations and measurement 
operations, we can always apply all the unitary operations first, and then apply 
the measurements, by the "Principle of Deferred Measurement" Chapter 4]. 
Then the unitary operations compose to a single unitary operation, and the mea- 
surements combine to form a single measurement operation Exercise 2.57]. A 
unitary operation U followed by a measurement operator {M„i} is the measurement 
operator {UMm}^ which also cannot perform cither transformation. So there is no 
physical process P to perform either transform. □ 

4.3. Quantum convolution and correlation are impossible. Finally, we can 
prove the main theorems and p^ : 

Proof of theorem Suppose such a process P existed. Then the composition 
{IQFT (g) IQ FT) o P o Q FT is physically realizable, which performs transformation 
p6| in theorem a contradiction. 

Thus there can be no physically realizable quantum convolution process. □ 

Note: this proof could use the fact that P would be linear, as would the compo- 
sition, which would contradict lemma ^ 

Proof of theorem The proof is the same as the proof of theorem |l^ with minor 
changes. Thus there can be no physically realizable quantum correlation process. 

□ 
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5. Conclusion and open problems 

We have shown there is no physical way to compute the convolution or correlation 
on quantum states. This is unfortunate for many applications, yet there may remain 
ways to attack specific sequences which yield efficient quantum algorithms. For 
example, it is possible to use a measurement to compute the transform ^ aibj\ij) — > 
X'^ aibi\ii) , but the probability of obtaining this state tends to as iV — > oo, likely 
rendering it useless in algorithms. 

A final interesting point is that given a sequence of numbers on a classical com- 
puter, we could compute a convolution, and prepare the quantum states in lemma 
p^ , which seems to violate theorem The point is that once the input sequences 
are in the quantum states, without knowing the classical information so we can 
make copies, there is no way to get the output state. Quantum states cannot be 
copied by the No Cloning Theorem |l^ , otherwise the transform of theorem |l^ 
could be computed by sampling enough copies of the input to determine the classi- 
cal sequences, classically computing the convolution, then constructing the output 
state. This approach is prevented by the No Cloning Theorem. The key point is 
that storing a classical information theory sequence of numbers into the quantum 
state coefficients is not a reversible process, since there is no way to read all the 
coefficients back out of the state. 

A final note: this result was inspired by a comment made by David Meyer, who 
obtained similar results independently. 
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